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Recently, Majorana Fermions (MFs) have attracted intensive attention due to their exotic statis¬ 
tics and possible applications in topological quantum computation (TQC). They are proposed to 
exist in various two-dimensional (2D) topological systems, such as px + ipy topological supercon¬ 
ductor and nanowire-superconducting hybridization system. In this paper, two types of Majorana 
Fermions with different polygon sign rules are pointed out. A “smoking gun” numerical evidence to 
identify MF’s classification is presented through looking for the signature of a first order topological 
quantum phase transition. By using it, several 2D topological superconductors are studied. 


I. INTRODUCTION 

Majorana fermion is a real fermion that is its own 
antiparticle[l|-Q. Because of its exotic properties and 
the possible exotic statistics j 3 -[T^. in condensed matter 
physics, the search for Majorana fermions (MFs) has at¬ 
tracted increasing interests. A variety of schemes to real¬ 
ize MFs (more accurately, Majorana bound states) have 
been proposed. A possible approach is to create a quan¬ 
tized vortex (tt-Aux) in the px + ip^-wave topological su¬ 
perconductors (SC) that traps MFs in vortex-core[2J-[^. 
Then, it is known that the quantized vortex in two dimen¬ 
sional (2D) topological SC with nonzero Chern number 
hosts a MF with exact zero energy. Another different 
approach to realize MFs is to consider a one-dimensional 
(ID) electronic nano-structures proximity-coupled to a 
bulk superconductor 0 , of which the unpaired Majo¬ 
rana fermions appear as the end-states. Then, based on 
this idea, several schemes are proposed to realize MFs 
that appear as the end-states of line-defects in 2 D non- 
topological SCs[ 29 |. 

To describe the Majorana zero mode, a real fermion 
field called Majorana fermion 7 = / dr[i^o' 0 * + 

(q/t = q/) is introduced. We consider a 2D gapped SC 
with a pair of Majorana modes with nearly zero ener¬ 
gies, of which the corresponding MFs are denoted by 
(715 72)- To describe the subspace of the system with 
two nearly degenerate states, the Fermion-parity opera¬ 
tor P = —27172 is introduced. Since = 1, P has two 
eigenvalues ± 1 , called even and odd Fermion-parities, re¬ 
spectively. Generally, there exists the coupling between 
two MFs and the effective Hamiltonian is given by 2^7172 
where t is the coupling constant. The quantum systems 
with multi-MFs (we call this lattice model to be Majo¬ 
rana lattice model) show nontrivial topological proper¬ 
ties, including a nonvanishing Chern number, chiral Ma¬ 
jorana edge statejsol-ls^. It was pointed out that the 
Majorana lattice model is really an induced ’’topological 
superconductor” on the parent TSC. 

A question arises, ’’Do MFs in different topological su¬ 
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perconductors belong to the same class?^^ Without no 
detectable degree of freedoms (MFs have zero energy, 
zero charge, zero spin), it is believed that MFs have ’’no 
HaiF\ Therefore, it was believed that all MFs in differ¬ 
ent models are same and belong to the same class. In this 
paper we point out that in 2D topological systems, MFs 
do have ” hair”, that is the number of quantized vortex, 
a topological degree of freedom. As a result, there exist 
two universal classes of MFs: MFs binding a tt-Aux and 
those with no flux-binding. We then introduce a topo¬ 
logical value to characterize the two classes of MFs and 
propose an (indirect) numerical approach to identify the 
class of MFs in topological systems. Our basis is the fact 
that multi-MFs with tt-Aux may be topologically different 
from those with no flux-binding in the Fermion-parity of 
the ground states and the order of quantum phase transi¬ 
tions. By calculating a dimensionless parameter, the gap 
ratio (see discussion below), the quantum number of the 
MFs becomes an observable ’’quantity”. 


II. CLASSIFICATION OF MAJORANA 
FERMIONS 

We begin by giving the definitions of two classes of 
MFs in 2D topological systems. One class contains the 
usual MFs without flux-binding (we call this class normal 
MFs); the other class contains composite objects of an 
MF together with a tt-Aux (we call this class topological 
MFs). See the illustrations of the two classes of MFs in 
FiglTl(a) and Figl 3 l(a). 

We consider a system with coupled MFs. The coupling 
strength is just the energy splitting from the intervortex 
quantum tunneling. We call this effective description as 
the Majorana lattice model, of which the tight-binding 
Hamiltonian can be written as 

P^'Tn.f — ^ ^ ^jktjkUkUj (^) 

(iA) 

where tjk is the hopping amplitude from j to k, and satis¬ 
fies tjk = tjj^. 7j is a Majorana operator (7! = jj) obey¬ 
ing anti-commutate relation {7j,7/c} = Sij = —Sji 

is a gauge factor. Thus, the total number of Majorana 
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modes N must be even and then we can divide the Majo- 
rana lattice into two sublattices. The pair (j, k) denotes 
the summation that runs over all the nearest neighbor 
(NN) pairs (with hopping amplitude t) and all the next- 
nearest neighbor (NNN) pairs (with hopping amplitude 
t'). Each triangular plaquette possesses ±7r/2 quantum 
flux effectively. 

This Hamiltonian allows different Z 2 gauge choices 
(sign rules) Sjk = ±1. From the theory of projective 
symmetry group, there are two possible gauge choices 
(sign rules) that correspond to two different classes of 
MFs: topological MFs obey a topological polygon sign 
rule; normal MFs obey a normal polygon sign rule. Ac¬ 
cording to topological polygon sign rule for topological 
MFs, there exists an extra phase related to a closed path 
that forms a polygon, given by half the sum of the in¬ 
terior angles of the polygon (ssf. Thus, for 4A^ (4A^ ± 1) 
topological MFs on a ring, there exists Ttt (±7r/2) flux 
inside the ring; for 4N ± 2 topological MFs on a ring, 
there is no flux inside the ring. On the contrary, accord¬ 
ing to normal polygon sign rule for normal MFs, there 
also exists an extra phase related to a closed path that 
forms a polygon. However, for AN ± 1 normal MFs on a 
ring, there exists ±7r/2 flux inside the ring and for AN or 
AN ± 2 normal MFs on a ring, there is no flux inside the 
ring. 


III. QUANTUM PHASE TRANSITIONS IN 
DIMERIZED MAJORANA RINGS 

A. The Hamiltonian of Majorana rings 

We begin from a (dimerized) Majorana ring of AN 
MFs (N is a positive integer number). The Hamiltonian 
reads 0 

2N — Si^N 

'HmF = i ^ {Ji7jb7j+l,a + Jlljaljb ( 2 ) 

i=i 

+ ^27i67j+l,6 + J2ljalj+l,a) 

where a, b denote the sublattices in a unit cell, Ji is the 
coupling constants between two MFs in a unit cell, J[ 
(J 2 ) are the (next) nearest coupling constants between 
two MFs in different unit cells. For the case of = 1, 
^i,Ar = 1; For the case of > 1, 6i^n = 0. For usual 
systems, we have J 2 < Ti, Tj. 

For a Majorana ring with AN normal MFs (we call it 
normal Majorana ring), due to normal polygon sign rule, 
we always have periodic boundary condition; while for a 
Majorana ring with AN topological MFs (we call it topo¬ 
logical Majorana ring), due to topological polygon sign 
rule, we have anti-periodic boundary condition owing to 
an extra tt flux inside the ring. FiglTJ(a) and Figl3l(a) 
illustrate a Majorana ring with 8 normal MFs and that 
with 8 topological MFs, respectively. 


B. Quantum phase transitions in normal Majorana 

rings 



FIG. 1: (color online) (a) Illustration of an normal Majorana 
ring with 8 normal Majorana Fermions (red dots); (b) Phase 
diagram of normal Majorana ring in thermodynamic limit, 
N 00 . The blue line denotes the first order topological 
quantum phase transition that switches Fermion-parity of the 
ground state and the red lines denote the second order phase 
transitions. 

Firstly, we study a (dimerized) Majorana ring with AN 
normal MFs. We pair {jjai Ijb) into a complex fermion 
asTja = (ci+cj)/\/2, 7j6 = -i{cj -j^/2, where Cj (c]) 
annihilates (creates) a complex fermion. Then the Ma¬ 
jorana ring’s energy spectra can be obtained through a 
fourier transformation Ck = In ther¬ 

modynamic limit, N ^ oc, the Hamiltonian in momen¬ 
tum space takes the form of 

Hmf = cos/c + Ji)N (3) 

k 

+ {J[ry - 2J 2 ) sin k]^k 

where = (c^,c_/c) and N, are Pauli matrices. The 
energy spectra are given by 

E±{k) = — 2 J 2 sink ± (J{)^ — 2Ji cos k Jf (4) 

where k = irn/N, n = 1,2,..., 2N. 

To characterize the topological properties of the 
ground states, we define a Fermion-parity operator in 
the following form 

2N 

p = (5) 

i=i 

Since P^ = 1, P has two eigenvalues ±1, called even and 
odd Fermion-parities, respectively. Due to [P, Hmf] = 0, 
the ground state \vac) should have a determinant parity 
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FIG. 2: (color online) (a) Illustration of the normal Majorana 
ring with four Majorana Fermions; (b) Phase diagram of the 
normal Majorana ring with four Majorana Fermions. The 
blue line denotes the first order topological quantum phase 
transition that switches Fermion-parity of the ground state; 

(c) The gap-ratio of the normal Majorana ring with four Ma¬ 
jorana Fermions. The divergence of the gap-ratio, a ^ oo, 
denotes a first order topological quantum phase transition; 

(d) The energy levels of the normal Majorana ring with four 
Majorana Fermions at J 2 = 0.4J{. The level-crossing denotes 
a first order topological quantum phase transition. 


P \vac) = ± \vac). For the case of P \vac) = — \vac ), 
the ground state is a topological phase; for the case of 
P \ vac) = \vac) , the ground state is a non-topological 
phase. 

For the Majorana ring described by Eq.([2]), the eigen¬ 
values of P is equal to 

sgn{E{k = 0) • E{k = tt)) (6) 


fers from the case in thermodynamic limit (or N oc) 
(See FiglTl(b)). The blue line in Fig 121 (b) denotes the 
first order TQPT that switches the fermion-parity of the 
ground state. To make the TQPT in a Majorana ring 
with only 4 MFs more clear, we plot the energy levels for 
the case J 2 = 0.4J{ in Figl2l(d). Now, the four energy 
levels are E±(k = 0), E±(k = tt). The level-crossing in 
Figl2l(d) shows the first order TQPT corresponding to 
the blue line in Figl2l(b)). 

In addition, to characterize the TQPT, we introduce a 
dimensionless parameter - gap-ratio^ 

_ maxl^fel - minl^fel 
minlE/cl 


where mdix\Ek\ and min|£’/c| are the maximum value 
and minimum value of the energy levels of the multi- 
MFs, respectively. When the gap-ratio turns to infinite, 
the energy gap closes and TQPT occurs. In the ther¬ 
modynamic limit, the gap-ratio a becomes where 


Wmf = maxlE/cl — min|£’/c| is the band-width and 
Amf = min \Ek\ is the energy gap of the Majorana ring. 
This is why we call a gap-ratio. On the other hand, for 
an normal Majorana ring with four MFs (or A = 1), due 
to max|£’/c| = \E{k = 7r)\ and min|E'/c| = \E{k = 0)\, 


the gap-ratio a is equal to 


E(k = 7T) 

Elk=0) 


— 1. In particular. 


at TQPT, accompanied by gap-closing, Amf = 0 or 
E(k = 0) = 0, the gap-ratio diverges. 


ex —y 00 . (8) 

See the results in Figl2l(c). 


C. Quantum phase transitions in topological 
Majorana rings 


where E{k = 0) = Ji — at /c = 0 and E{k = tt) = 
JiE Ji k = E. Thus, at Ji = Ji, the energy gap closes, 
E±{k) = 0, at which a topological quantum phase tran¬ 
sition (TQPT) occurs 0. In Fig|TJ(b), we plot the global 
phase diagram of the normal Majorana ring in thermo¬ 
dynamic limit, N ^ 00 . There are three phases: gapless 
phase (we don’t discuss this phase due to triviality), topo¬ 
logical phase and non-topological phase. The blue line 
denotes the first order TQPT that switches the Fermion- 
parity of the ground state 0 and the red lines denote the 
second order phase transitions. In yellow region (Ji > 

J 2 < Ji)i we have P = 1 and the ground state cor¬ 
responds to a trivial state (non-topological phase) with 
even Fermion-parity, P \vac) = \vac). In green region 
(Ji < Ji, J 2 < Jj), we have P = — 1 and the ground state 
corresponds to a topological phase with odd Fermion- 
parity, P \vac) = — |uac)0. So, for a normal Majorana 
ring with arbitrary A, a first order TQPT that switches 
the Fermion-parity occurs at Ji = J^. 

On the other hand, for a Majorana ring with 4 normal 
MFs (or N = 1), the phase diagram (See Figl2l(b)) dif- 


Next, we study a (dimerized) Majorana ring of 4A 
topological MFs. Due to the extra tt-Hux inside the ring, 
we have an anti-periodic boundary condition of the topo¬ 
logical Majorana ring. The energy spectra are the same 
to those of the normal Majorana ring as 

E±(k) = - 2 J 2 sink ± f {J{r - 2J{Jicosk + Jf (9) 

with different wave-vectors, k = 7r(n — l/2)/A, n = 
1 , 2 ,..., 2A. Now, there are no high symmetry points at 
k = 0^ k = E. As a result, the ground state always has 
even fermion-parity, P \ vac) = \vac). In thermodynamic 
limit, a second order phase transition occurs at Ji = Ji 
between a non-topological phase (the left yellow region in 
Figl3l(b)) and another non-topological phase (the right 
yellow region in Figl3l(b)). For the topological Majorana 
rings with finite A, the energy levels can smoothly change 
from one phase to the other without level-crossing. At 
the quantum critical point, Ji = J^, in the thermody¬ 
namic limit, we have a ^ N ^ oo; while for finite A, a 
is a finite value. 
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logical Majorana ring with four MFs are E±{k = 7 r/ 2 ), 
E±{k = — 7 r/ 2 ). In FiglU(d), the energy levels of the 
topological Majorana ring with four Majorana Fermions 
are plotted for the case of J 2 = O.lJj. The gap-ratio for 
topological MFs is obtained to be 


E-^{k = — 7 r/ 2 ) 
E_^{k = 7r/2) 


For the case of J 2 ^ Ti, Tj, due to E±{k = 7 r/ 2 ) 
E±{k = — 7 r/ 2 ), we always have a very small gap-ratio as 


a ^ 0. 


( 10 ) 


FIG. 3: (color online) Illustration of a topological Majorana 
ring with 8 topological Majorana Fermions (red dots with 
strings). Due to the polygon rule, there exists an extra tt-Hux 
inside the ring for topological Fermions; (b) Phase diagram of 
topological Majorana ring in thermodynamic limit, N 00 . 
The red lines denote the second order phase transitions. 



FIG. 4: (color online) (a) Illustration of the topological Majo¬ 
rana ring with four Majorana Fermions. Due to the polygon 
rule, there exists an extra tt-Hux inside the ring for topologi¬ 
cal Majorana Fermions; (b) Phase diagram of the topological 
Majorana ring with four Majorana Fermions. The blue line 
denotes the first order topological quantum phase transition 
that switches Fermion-parity of the ground state; (c) The gap- 
ratio of the topological Majorana ring with four Majorana 
Fermions; (d) The energy levels of the topological Majorana 
ring with four Majorana Fermions at J 2 = 0.1J{. 


On the other hand, for a Majorana ring with 4 topo¬ 
logical MFs (or N = 1)^ the phase diagram differs from 
that for the normal case. See the results in Figlll(b). 
A first order TQPT occurs in the region with large J 2 
that is irrelevant to traditional topological systems. In 
the region with small J 2 , the first order TQPT will never 
occur. In Figlll(c), the four energy levels for a topo- 


IV. NUMERICAL METHOD TO IDENTIFY 
THE MAJORANA FERMION’S 
CLASSIFICATION 

From above discussion, the sharp distinctions between 
normal and topological MFs are found in relevant physics 
(J 2 < Ji, J 2 < T{): for four normal MFs on a ring, at 
the point of first order TQPT, a ^ 00 (or a ^ 1); for 
four topological MFs on a ring, without the first order 
TQPT, (a ^ 0 (or a < 1). Therefore, we propose a 
numerical method to distinguish the normal/non-normal 
polygon sign rule for the MFs by calculating the gap- 
ratio (a in a topological system with four MFs that form 
a (dimerized) Majorana ring. 



FIG. 5: (color online) Illustration of the varying of a Majorana 
ring with four Majorana Fermions from the limit LxfLy ^ 1 
(or Ji/Ji > 1) to the limit LxjLy > 1 (or Ji/Ji 1). 

In the first step, we study the given 2D topological 
system without considering the MFs. After diagonaliz¬ 
ing the BdG equation, we obtain the energy spectra in 
momentum space F^(k) and the energy gap of the system 

In the second step, the energy levels of the 2D topo¬ 
logical system with a pair of MFs ( 71 , 72 ) are calculated 
by numerical approach. We can derive the energy lev¬ 
els of the MFs with almost zero energies, EE^ {E > 0). 
When there are two MFs nearby, the quantum tunneling 
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effect occurs and the two MFs couple. The energy split¬ 
ting between two nearly zero modes /S.E = 2 E versus the 
distance L of the two MFs can be obtained. In general, 
AE" oscillates and decreases exponentially with L and can 
be described by a function as AE oc |cos(E • /cf)| 

where ^ ^ vpl^f is the correlated length ^ and vp is 
the Fermi velocity. We plot the enveloping line of AE 
(we denote it by AEei) by choosing the distance to be 
Eel ^ lYn/kp where n is an integer number. It is obvi¬ 
ous that AEei becomes a monotonous function via L and 
decays exponentially, AEei oc . 

In the third step, we consider the topological system 
with four localized MFs that form an x Ly square (a 
dimerized Majorana ring). The distance between 71^, 715 
and that between 72a, 726 is E^,, the distance between 71^, 
726 and that between 715, 72a is Ly. If we fix Ly (or fix 

), and then we can smoothly tune J\ by changing L^ = 
Eel ^ 'Knjkp. In the limit of E^, ^ 00, we may have J\ 
in the limit of E^, ^ 0, we may have Ji':^ J^. See the 
illustration in FiglSl During varying L^jLy (or Ei/Ei and 
eventually a), we carefully look for the signature of a first 
order TQPT with level-crossing that indicates a diverge 
gap-ratio, or a ^ oc. For a 2 D topological system on 
lattice, due to the discreteness of L^jLy (or J\j the 
gap-ratio a will never diverge but may be a fairly large 
value. The large gap-ratio (for example, max a > 10 ) 
could be regarded as a strong evidence of the first order 
TQPT. Eventually, the MFs obey normal polygon sign 
rule. On the contrary, if the resulting gap-ratio a in a 2D 
topological system is always a small value (for example 
max (a < 1), we can exclude the possibility of a first order 
TQPT and ultimately identify the topological polygon 
sign rule of the MFs. 


V. CLASSIFY MAJORANA BOUND STATES IN 
2D TOPOLOGICAL SUPERCONDUCTORS 

In condensed matter systems, MFs are proposed to 
exist in various two-dimensional (2D) topological SCs [ 3 - 
flsl . In 2D strong topological SCs (so termed because of 
the non-zero Chern number), MFs could be induced by 
quantized vortices or dislocations [ 35 [ [ 36 |; in 2 D weak 
topological SCs (so termed because the Chern number is 
zero), MFs could be induced by line defects [ 3 , [ 29 |. In 
the following, we studied MFs in strong topological SCs 
in Sec.III.A and B and those in weak topological SCs in 
Sec.III.C and D. 


A. Majorana bound states in a 2D px + ipy 
topological superconductor 

In the first example, we studied the MFs around vor¬ 
tices in d. 2 D Px E ipy topological SC. The Hamiltonian 



FIG. 6: (color online) (a) Illustration of two Majorana bound 
states trapped by 7r-fluxes in px + ipy topological supercon¬ 
ductor; (b) The particle density distribution of Majorana zero 
modes around two 7r-fluxes; (c) The energy splitting AEei via 
E (the distance between the two Majorana bound states); (d) 
The gap-ratio. 


of a Pa, + ipy SC on a square lattice is written as[ 25 | 
Hpip = E 

j tl=x,y j 

j 

where cj is an electronic annihilation operator, fi is the 
chemical potential, A is the SC pairing-order parameter 
and t is the hopping strength. The lattice constant was 
set to unity in this paper. In the following, we chose the 
parameters ast = l,/i = —1,A = 0 . 4 . The ground state 
was the weak-pairing phase that is a (strong) topological 
SC. 

We first studied two MFs (71, 72) around two vortices 
with numerical calculations on a 120 x 40 lattice. We 
found that there exists a Majorana zero mode around 
each vortex. Note the particle density around two tt- 
fluxes in Figl6l(b). When there are two fluxes nearby, 
inter-flux quantum tunneling occurs and the two MFs 
couple. Fig El (c) shows the energy splitting AEei. Next, 
we studied four coupled MFs of vortices, 71^, 715, 72a, 
726, that formed an Lx x Ly square (a dimerized Ma¬ 
jorana ring). See inset in FigEl(d). By fixing Ly at 6 
and varying E^,, we calculated the gap ratio a and show 
the result in FigEl(d), in which one can see a very small 
gap ratio. There is no evidence of the first order TQPT. 
Therefore, we identified the MFs induced by the vortices 
inpxEipy topological SCs to be topological MFs. It is ob¬ 
vious that this conclusion (topological MFs in 2D p^ipy 
topological SCs) is consistent with earlier results!^. 
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B. Majorana bound states in an s-wave topological 
superconductor with Rashba spin-orbital coupling 



(a) (b) 


FIG. 7: (color online) (a) The energy splitting AF^ei via the 
length L of two vortices. The inset in (a) shows the particle 
density distribution of the Majorana bound states; (b) The 
gap-ratio of four Majorana bound states. The inset in (b) 
shows two positive energy levels of four coupled Majorana 
bound states. 

The second model is MFs in an s-wave SC with Rashba 
spin-orbital (SO) coupling on a square lattice [6|-[8|, [s^- 
The Hamiltonian is given by i7s-wave-SO = i^kin + ^so + 
Hsc where the kinetic energy term F^kin, the Rashba SO 
coupling term Hsoi and the SC pairing term are given 
as 

^kin = (^j+Fcr^icr + cj_^^Cjcr) (12) 

jo- ^=x,y 

- - cl^Cji), 

jo- j 

Hso = 

j 

+ + H.C., 

4c = -A^(cj,4 + if.c.). 

j 

Here, Cja- annihilates (creates) a fermion at site 

j = {jx, jy) with spin a = (t, I), /i = x or g', which 
is a basic vector for the square lattice. A serves as the 
SO coupling constant and A as the s-wave SC pairing- 
order parameter, fi is the chemical potential and h is the 
strength of the Zeeman field. The lattice constant was 
set as unity. The parameters were chosen to be t = 1, 
fi = —4, h = 0.8, A = 0.5, A = 0.5. In this case, the 
ground state was a (strong) topological SC. 

We then studied two MFs around two vortices through 
numerical calculations on a 100 x 24 lattice. The parti¬ 
cle distribution of the Majorana zero mode is given in 


the inset in Figl3(a). The results of the energy splitting 
AE'ei via the distance of the two vortices L are given in 
Fig 13 (a). Next we studied four coupled MFs around the 
vortices 71 ^, 715 , 720 , and 725 , that formed an x Ly 
square (a dimerized Majorana ring). We fixed Ly to be 
6 and varied Lx. The gap ratio is shown in FigI71(b). 
One can see a small gap ratio. As a result, we con¬ 
clude that MFs in an s-wave topological superconductor 
with Rashba SO coupling obey topological polygon sign 
rulellij. 

C. Majorana bound states in a nanowire-SC 
hybridization system 

The third model was a ID semiconducting nanowire 
with strong spin-orbital coupling in a Zeeman field, 
proximity-coupled to an s-wave superconductor. See the 
illustration in Figl 8 l(a). The Hamiltonian of the system 
is ^nano-SC = ^iD + ^sc, where ^ID describes a ID 
semiconducting nanowire; it is written as 

^ID = —ts ^ + C^_^^Cjcr) + H.C. (13) 

jcr 'fi=x 

- 

j j 

3 jo 

Hsc describes the 2D superconductor out of ID semicon¬ 
ducting nanowire and is written as 

H.C = -tJ2Yl 

jo y=x,y 

- M T! - A 

jo j 

where Cjcj is an electronic annihilation operator and t (t^), 
yU, A, A, h denote the hopping parameters, the chemical 
potential, the (induced) pairing order parameters, the 
spin-orbit coupling strength and the Zeeman field, re¬ 
spectively. The lattice constant was set to be unity. In 
the following, we chose the parameters as G = 1 , = 0 . 8 , 

fi = — 2 , A = 0.5, A = 0.4, t = 0.4. The ground state of 
the system is a (weak) topological SC. 

Then, treating a nanowire by numerical calculations on 
a 120 X 20 lattice, we found that there are two Majorana 
zero modes ( 71 , 72 ) at two ends of a nanowire. See the il¬ 
lustration in Figl 8 l(a). There exist two zero modes in the 
energy gap shown in Figl 8 l(b). In Figl 8 l(c), we plot the 
particle distribution of the zero modes. The relationship 
between the energy splitting A£)ei and the length L of 
the nanowire is shown in Figl 8 l(d). Next, we considered 
two parallel ID semiconducting nanowires on an SC. See 
the right inset in FiglQl Here the low energy physics is 
dominated by four coupled MFs, 71 ^, 715 , 72 a, 726 , that 
form an Lx x Ly square (or a dimerized Majorana ring). 
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FIG. 8: (color online) (a) Illustration of nanowire-SC hy¬ 
bridization system; (b) Two Majorana bound states with zero 
energy; (c) Particle distribution of Majorana bound states at 
the ends of the nanowire on s-wave SC; (d) The energy split¬ 
ting of two Majorana bound states. The inset in (d) shows a 
very tiny J 2 . 



FIG. 9: (color online) The gap-ratio of the four Majorana 
bound states in nanowire-SC hybridization system. The right 
inset illustrates two semiconducting nanowires on an SC sys¬ 
tem. The left inset shows two positive energy levels of four 
Majorana bound states. 


We fixed the distance between two parallel nanowires to 
he Ly = 2 and then changed the length of the two 
nanowires. The left inset of Figl9] shows the two positive 
energy levels from four MFs of the two nanowires vs. 

In Figini the gap ratio a is obtained. From FiglHl one can 
see that the maximum value of a reached 17 at = 30. 
The sharp enhancement of the gap ratio a is obviously 
the consequence of a first order TQPT. As a result, we 
conclude that the MFs in the nanowire-SC hybridization 
system obey normal polygon sign rule. 


D. Majorana bound states in a p-wave 
superconductor on a honeycomb lattice 

The fourth model was MFs in a 2D p-wave supercon¬ 
ductor on a honeycomb lattice. The Hamiltonian of a 
p-wave superconductor for spinless fermions on a honey¬ 
comb lattice is written as|^ 

-^p—wave ~ t ^ ^ Cj t ^ ^ C ^ ^ H.C.^ 

(*i> ((*i)> (*i) 

(15) 

where t (t') denote the strengths of nearest (next nearest) 
neighbor hopping. The p-wave pairing order parameters 
are defined by = A, A^j+as = 0- 

^oi (<^ = 152 , 3 ) denotes a vector that connects the near¬ 
est neighbor sites i and i + Along the red links in 
FigllQl(a) and Fig(T0l(c), the SC order parameter is fi¬ 
nite; along black links, the SC order parameter is zero. 
The lattice constant is set to be unity. We chose t = 1, 
t' = 0.01, A = 1.34 in this section. Here the ground state 
is a (weak) topological SC. 



FIG. 10: (color online) (a) Illustration of two Majorana bound 
states induced by a line defect; (b) The energy splitting AF^ei 
via the length L of the line defect; (c) Illustration of four 
Majorana bound states induced by two line defects; (d) The 
gap-ratio of the four coupled Majorana bound states induced 
by two parallel line defects. 

In Ref. [ 23 , it was found that there exist two Majorana 
modes localized at the ends of the line defect (LD), 71 , 72 . 
See the illustration in FigJTOKa). The end of an LD can 
be considered to be the boundary of a one-dimensional 
p-wave SCQ. Thus, each end of the LD traps a dangling 
Majorana zero mode. We studied two MFs around a LD 
with numerical calculations on a 100 x 50 lattice and give 
the results of the energy splitting AE’ei via the length L 
of the LD in Fig(T0l(b). Then, as shown in FigJTOKc). 
we studied four coupled MFs of two parallel LDs that 
formed an x Ly square (a dimerized Majorana ring). 
The distance between the two parallel LDs was fixed to be 
2 (or Ly = 2). We varied the length of the LDs, L^. The 
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results of the gap ratio a are given in Fig(T0l(d), in which 
the maximum value of a reaches 6.5 at = 9. These 
results indicate a first order TQPT with level-crossing. 
Thus, the MFs in p-wave superconductors on honeycomb 
lattice also obey normal polygon sign rule. 

VI. DISCUSSION AND CONCLUSION 



Normal Majorana Fermions 

Topological Majorana Fermions 

Definition 

Pure Majorana ^ 
Fermions 

Maj orana F ermions ^ ^ 

with a n phase string 

Statistics 

Non-Abelian 

Non-Abelian 

Polygon 

rule 

Normal polygon sign rule 

Topological polygon sign rule 

Examples 

Line-defect-induced MBSs 
in 2D weak topological SCs: 

1) nanowire-SC 
hybridization system; 

2) p-wave SC on honeycomb 
lattice. 

Vortex-induced MBSs in 2D 
strong topological SCs: 

1) px+ipy topological SC; 

2) s-wave topological SC with 
Rashba spin-orbital coupling. 


FIG. 11: (color online) The difference between normal Majo- 
rana Fermions and topological Majorana Fermions. 


In the end, we draw our conclusions (see the sum¬ 
mary in Fig. 11). We have pointed out that in 2D 
topological superconductors, there exist two classes of 
MFs: MFs obeying normal polygon sign rule and MFs 
obeying topological polygon sign rule. A numerical ap¬ 
proach was proposed to identify the polygon sign rule 
of the MFs by looking for the signature of a first order 
TQPT of multi-MFs. Applying the approach to study 
several 2D topological systems, we found that vortex- 
induced MFs in 2D strong topological SCs (a Px + ipy 
topological superconductor in Sec.VI.A and an s-wave 
topological superconductor with Rashba spin-orbital cou¬ 
pling in Sec.VLB) obey topological polygon sign rule 
and line-defect-induced MFs in 2D weak topological SCs 
(nanowire-SC hybridization system in Sec.VI.C, p-wave 
superconductor on honeycomb lattice in Sec.VI.D) obey 
normal polygon sign rule. 
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